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Abstract 

We introduce a special class of truncated Weyl-Heisenberg algebra and discuss the corresponding 
Hilbertian and analytical representations. Subsequently, we study the effect of a quantum network 
of beam splitting on coherent states of this nonlinear class of harmonic oscillators. We particularly 
focus on quantum networks involving one and two beam splitters and examine the degree of bipartite 
as well as tripartite entanglement using the linear entropy. 



* daoud@pks.mpg.de, m_daoud@hotmail.com 
^ ajellal@ictp.it, jellal@ucd.ma 
'choubabi@gmail.com 
^ elkinani@gmail. com 



1 Introduction 



Over the past few years, the entanglement of quantum systems was reahzed to be a valuable and 
crucial resource in quantum information processing. It allows for powerful new communication and 
computational tasks that are not possible classically. One may quote for instance quantum teleporta- 
tion [T], superdense coding [2], quantum key distribution [3], telecloning [3], quantum cryptography [S] 
and quantum computation [6l [7] . Then it is not surprising then that over the last decade much efforts 
has been devoted to capture, quantify and assess the power of quantum entanglement. In this sense, 
many authors studied the development of a quantitative theory of entanglement and the definition of 
its basic measure (concurrence, entanglement of formation and linear entropy [Sll^limillj). Entangled 
quantum systems can exhibit correlations that cannot be explained on the basis of classical laws and 
the entanglement in a collection of states is clearly a signature of non-classicality |12j . 

In quantum theory, the states minimizing the quantum fluctuations, which are closest to the classi- 
cal ones are the coherent states [l'S\ [Hj . This motivated the considerable interest in the entanglement 
of coherent states \15 \ [TT } I18j. The coherent state approach is not just a mathematical tool, but it 
my also helps to understand the entanglement for such states, which provide a bridge between quan- 
tum and classical worlds. It is also important to stress that in connection with quantum entanglement, 
many experimental results were obtained. In fact, the experimental generation and characterization of 
the entangled electromagnetic states can be achieved using type I or type II parametric down conver- 
sion [19]. Another experimentally accessible device, which can be used to generate optical entangled 
states, is the beam splitter [201 ED [22] ■ In quantum optics the action of a beam splitter, which is 
essentially a coupling of two electromagnetic modes, can be represented by a unitary operator relating 
the input and the output states. In general, the output state is a superposition of the Fock states 
which is entangled, except the harmonic oscillator coherent states. Indeed, the harmonic oscillator 
coherent states does not exhibit entanglement when passed through one arm of 50:50 beam splitter 
while the second arm is left in the empty vacuum state [23]. In the same sense, the entanglement 
behavior of the SU{2) spin coherent states, when passed through a beam splitter, was examined in 
[12j . Similar study was done in [24J for the SU{1, 1) coherent states. 

On the other hand, the truncated harmonic oscillator was used by Pegg and Barnett [25] to 
define the phase states for the quantized single modes of the electromagnetic field. They suggested 
to truncate to some finite (but arbitrarily large) order the infinite-dimensional representation space 
of the oscillator algebra. This was done to get rid of the difficulty related to the infinite-dimensional 
character of the representation space of the Weyl-Heisenberg algebra, which constitutes a drawback 
in defining a phase operator in a consistent way |26t [271 128j . Motivated by these investigations and 
in particular [25] , we propose a refined version of the truncated oscillator algebra introduced by Pegg 
and Barnett. More precisely, we introduce a nonlinear class of Weyl-Heisenberg algebra and analyze 
its representations. We construct the associated coherent states using a suitable analytical realization 
and we study the degree of entanglement of such states when passed through a quantum network of 
beams splitters. 

The outline of the paper is as follows. In section 2, we introduce a generalized Weyl-Heisenberg 
algebra. We discuss the corresponding Hilbertian representation and the analytical Bargmann realiza- 
tion. It is remarkable that the obtained analytical realization provides us with an over-complete set of 
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states and then constitute a system of coherent states in the Klauder-Perelomov sense. In section 3, 
we define the action of a quantum network of k beam sphtters. We show that this action leads to the 
SU{k + 1) coherent states labeled by complex variables and related to the reflection and transmission 
parameters of beam splitters. Based on this result, we investigate in detail the effect of one and 
two beam splitters on the generalized Weyl-Heisenberg algebra coherent states. We derive the linear 
entropy to study the bipartite entanglement degree of the output states, study the entanglement in a 
system of three particles. Concluding remarks close this paper. 

2 Weyl— Heisenberg algebra and Fock— Bargmann realization 

We start by introducing and discussing some interesting properties of a generalized version of the 
truncated harmonic oscillator, which involves a half positive integer s. In the limiting case when this 
parameter goes to the infinity, one recovers the usual harmonic oscillator with infinite dimensional 
Fock space. This provides us with another mathematical tool to deal with truncation of the ordinary 
harmonic oscillator that is different from one discussed in |25) . We mention that the idea of truncated 
Weyl-Heisenberg is mainly inspired by the polynomial deformation of Lie algebras introduced in 
[50] and extensively discussed in the context of quantum algebras [311 [321 [33] . The truncated Weyl- 
Heisenberg algebra discussed here can be viewed as a special case of /-deformed oscillators introduced 
in \34:\ 135] which is relevant in the algebraic description of a large class of nonlinear quantum systems. 

2.1 Finite dimensional Weyl-Heisenberg algebra 

The basic ingredient that we shall use in what follows is the generalized Weyl-Heisenberg algebra 
characterized by a positive real parameter (s > 0). This algebra is generated by the set of operators 
{a'^,a^ ,N, 1} satisfying the relations 

[N,a-] = -a-, [N,a+] = +a+, [a",a+]=I-— (1) 

s 

where N is the number operator and the element I commutes with all other operators. Clearly, when 
s —7- c«, end up with the ordinary Weyl-Heisenberg algebra. In this respect, the parameter s can be 
regarded as a measure of the distortion of this algebra. For convenience, we assume that 2s E N. The 
generalized oscillator algebra can be naturally represented, on the Fock space T of the eigenstates of 
the number operator A^, by 

N \n) = n \n), {n\m) = 5nm, n, m G N (2) 

and the vacuum state satisfies a~|0) = 0. We define the actions of the operators a~ and a+ on J-" as 

a~\n) = y/F{n)\n-l), a+|n) = ^^((n + l)|n + 1). (3) 

Note that, a"*" and a~ are mutually adjoint, namely a+ = (a^)^ and N is, in general, different from 
the product a+a^. The structure function F(-) is an analytic function with the properties -F(O) = 
and F{n) > 0, n = 1, ■ ■ ■ . F{-) is characteristic to the distortion or the truncation scheme and satisfies 
the following recursion formula 

F{n + 1) - F{n) = 1 - - (4) 
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which gives by simple iteration the form 

71 

F(n) = — (2s + l-n). (5) 
2s 

Note that, the structure function obeys the condition 

F{2s + 1) = (6) 

and the creation-annihilation operators satisfy the nilpotency relations 

(a-)2^+i = 0, {a+f'+^ = 0. (7) 

This means that the corresponding representation is (2s + l)-dimensional. As evoked above the 
truncated oscillator algebra ([T]) constitutes a particular variant of the /-deformed oscillators [3H [35] . 
Indeed, one has 

a- = b-f{N), a+ = /(iV)6+, N = b+b- , (8) 
where the function f{N), reflecting the distortion from the usual bosons, is given by 



The equation ([8]) is a Holstein-Primakoff realization of the algebra ([1]). Finally, we point out that for 
s large, we have 

a± ~ b^ (9) 

traducing that the distorted algebra ([T|) coincides with the linear harmonic oscillator one. More 
importantly, the generalized Weyl-Heisenberg algebra provides us with a simply framework to deal 
with finite dimensional oscillator system. It is important to stress that this algebra is similar to one 
introduced in [3U] in order to define the phase operator for nonlinear oscillators and to derive the 
associated temporally stable phase states. In this respect, the generalized Weyl-Heisenberg algebra 
can be viewed as a refined version of truncated harmonic oscillator. 

Of course the problem of practical realization of the finite dimensional nonlinear quantum oscillator 
algebra arises. At this point one should emphasize that many exactly solvable systems enjoy the 
generalized Weyl-Heisenberg symmetry and are possessing finite dimensional representation space. 
One may quote for instance one dimensional quantum systems having finite discrete spectrums like 
ones evolving in the modified Poschl- Teller [37] and Morse [38] potentials, see also [39]. Quantum 
systems described by a nonlinear Hamiltonian are familiar in the context of nonlinear quantum optics 
as for instance electromagnetic field propagating trough a nonlinear Kerr medium [40] . Indeed, the 
dynamical evolution of the electromagnetic field propagating in a Kerr medium can be described by 
the following Hamiltonian 

HKerr = b+b' - K b+'^b'^ 

where k characterizes the Kerr nonlinearity. Setting ^ = the Hamiltonian Hxew coincides with the 
function structure F{N) with N = b^b^ . It is also interesting to stress that the ladder operators a"*" 
and a~ can be related to Stokes operators introduced in [HI |l2] in order to define a unitary operator 
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representing the exponential of the phase difference between two modes of the electromagnetic field. 
This relation can be expressed as 

s+ = y/sa+, S- = y/sa-, S3 = ^{N-s) 

and one can simply verify that the Stokes generators ,5+, S- and S3 satisfy the commutation relations 
of the su(2) algebra. 

We close this subsection by noting that, in recent years, a special interest is devoted to the pos- 
sibility of generating and manipulating systems, whose dynamics can be closed within a finite set of 
n-photon states. Such systems are very important for the implementation of models in the quantum 
information theory. The experimental devices leading to the finite-dimensional state generation are 
referred, in the literature, as linear |43t I44j or non-linear |45|, H6] quantum scissors. 

2.2 Fock-Bargmann realization 

The second ingredient needed for our task is the coherent states associated with the above generalized 
algebra. A simply way to construct these states is to use the analytical Bargmann representation. 
Indeed, we realize the annihilation operator 

a-^^ (10) 
az 

as derivation with respect to a complex variable z. The elements of the Fock space are realized as 
follows 

|n) CnZ^. (11) 

Using the action of the operator a~ on the Fock space F given by ([3]) and the correspondence (fTTj) . 
one obtains the following recursion formula 



(n + l)Cn+i = yjF{n+l)Cn. (12) 
It follows that the coefficients C„ are given by 



^" - ^°\/ (2s)"n!(2s-n)!- ^^^^ 

To simplify, we set hereafter Co = 1. Having the expression of the coefficients Cn, one can determine 
the differential action of the creation operator a"*". Indeed, since the operator N acts as 

N ^z^ (14) 
dz 

one can easily see, using the action of the generator on the Fock space, that 

2;2 d 

On the other hand, a general vector of F 

2s 



n=0 
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is represented in the Bargmann space as 

2s 



J]0nC„z". (16) 

ra=0 

Tlie inner product of two functions (j) and cp' is defined now as 



[z) = 

n=0 



d'z^{z)(l)'{z)(l){z) (17) 

where the integration is carried out on the whole complex plane. The integration measure S, assumed 
to be isotropic, can be computed by choosing \(p) = \n) and \(f)') = \n'). Thus, one has to look for a 
solution of the following moment equation 

f+°° , , , n„,i (2s)"n!(2s-n)! , , 

27r d£.S(£.)£,2«+i = ^—^-^^^ (18) 

where = |2:|. To find the function satisfying the equation (jlSp . we use Mellin transform technique to 
obtain 

E(,) = ^(l + ^)-2-^ (19) 
vr Zs 

At this stage, one can write the function (j){z) as the product of the state |0) with some ket \z) labeled 
by the complex conjugate of the variable z. This is 

cPiz)=Mm. (20) 

Taking {(p) = |n), we obtain 



2s 



\z)=My J-—^ -z^\n). (21) 

The normalization constant for the the states (j2ip is given by 

^•=(l + te)-*. (22) 

It is important to remark that the states (j2ip are coherent in the Klauder-Perelomov sense |131 [T3] ■ 
Indeed, it is easy to see that the states \z) can be written as 

\z) =7^exp(2a+)|0). 

They satisfy the over-completeness property 



/ cPzn\z\)\z){z\ = Y,\n){n\ 

71=0 



where the measure S(|2;|) is given by ()19p . It should be noticed the states ()2ip are similar to nonlinear 
coherent states derived in [Uj . They coincide with the standard Glauber coherent states for the usual 
harmonic oscillator at large s. 
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3 Quantum network of beam splitters 



In the recent years, the study of entangled states has revived interest in the beam sphtter. This is due 
to the fact that this device offers a simple way to probe the quantum nature of electromagnetic field 
through simple experiments. As mentioned before, many authors studied the behavior of quantum 
states when passed through a beam splitter. This device is an optical element with two input ports 
and two output ports that, in some sense, governs the interaction of two harmonic oscillators. The 
input and output boson operators are related by a unitary transformation which is an element of 
the SU{2) group. Recently, a quantum network of beam splitters was used to create multiparticle 
entangled states of continuous variables [IS] and also multiparticle entangled coherent states |49] . 
Here we examine the entanglement of the coherent states (|2ip passing through a quantum network of 
beam splitters. Note that, a set of k beam splitters can be experimentally used to generate SU {k + 1) 
coherent states labeled by the reflection and transmission parameters of the involved beam splitters. 

3.1 Generation of SU{k + 1) coherent states 

The most general transformation defining a quantum network of k beam splitters is given by the 
unitary operator 

Bk = I3k,k+i{ek)l3k-i,k{(^k-i) • • • i3i,2(^i) (23) 



where the operators 

'" i 
2 



Si,i+i{9i) = exp 



(24) 



are the ordinary SU{2) beam splitters with / = 1,2,- •• ,k. The objects 6^ and 6; are the usual 
harmonic oscillator ladder operators. The reflection and transmission coefficients 

ti = cos— , ri = sm—. (25) 

are defined in terms of the angles 9i. The operator Bk is actually acting on the states |ni, n2, • • • , Uk) 
of the usual /c-dimensional harmonic oscillator. If the input state is |ni, n2, • • • , Uk), then the Bk action 
leads to the following Fock states superposition 

Bk\ni,n2, ■■■ , Uk+i) = ^ iSlTi'C/'- ,nriY' \mi,m2, • • • , ruk+i) (26) 

and in general the output is a (A; + l)-particle entangled state. On the other hand, the action of the 
unitary operator Bk on the state |ni, 0, • • • ,0) gives 



^,ini,o,... ,0) = cyy--- y '"^'^ 

hi^o n.:r^oVK-2)!(n2-n3)!...(n,-n.,,)!n..,! 
x|ni - n2,n2 - 71-3, • • • ,nfc+i) (27) 



where the normalization constant reads as 

|2\-"i/2 



C = (1 + 161^ + 161^161^ • • • + ■ ■ ■ MT^'^' (28) 
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and the new variables ^ 

6 = itz+i^ for (/ = 1,2,--- & = (29) 

are defined in terms of the reflection and transmission coefficients of the beam splitters constituting 
the network. Then, the output state ()27p turns out to be the SU{k + l) coherent state associated with 
the completely symmetric representation labeled by the integer ni, see for instance [3U] where similar 
notations were used. This means that a network of k beam splitters can be used as experimental 
device to generate SU{k + 1) coherent states. 

3.2 Bipartite entanglement of truncated harmonic oscillator coherent states 

We start by investigating the effect of one beam splitter on the coherent states ([^T|) associated with 
the generalized or truncated harmonic oscillator introduced in the first section. As discussed above, 
for s large, the dimension of the generalized oscillator Fock space becomes infinite. In this case one can 
easily check that the coherent states ()2ip reduce to Glauber coherent states of the infinite dimensional 
harmonic oscillator. It has been proven that the Glauber states are the only pure states that when 
passed through one arm of the beam splitter, the output resultant states are disentangled [51]. In this 
respect, the main task here is to determine the behavior of the entanglement in terms of the Fock 
space dimension and to prove that the entanglement disappears when the Fock space dimension goes 
to infinity. For that end, we consider the situation where the coherent state ()2ip is injected into one 
port while no photon is injected into the other port. 

Using the equation (|27p for A: = 1, it is simply verified that the action of the beam splitter on the 
input state \z,0) can be cast in the following form 



2s 2s— n 



fi,.m)i..o) g ^ p,_i_'^^),„,,, (30) 

where M is given by (f22]) . To determine the degree of entanglement of the beam splitter output sates, 
we use the linear entropy which is the upper bound of the Von Neumann entropy [12]. This is 



S = l- Tt{pI) (31) 

where pi = Tr2/9i2 is the reduced density operator obtained by tracing over the states of the second 
system. Using ([30]) . the linear entropy is obtained as 

2s 2s min(2s— ?i,2s— n') min(2s— n,2s— n') 

"5 = 1 -EE E E h{n,n\q)h{n\n,ci) (32) 

n=On'=0 q=0 q'=0 

where 



b{n,n',q) = |AA|2(2s)!^^ ^n+n (33) 

^ ' '^^ ' '{23)^1 {./2^)n+n' \l {2s - n - qy.{2s - n' - q)\n\n'l ^ ' 

At this stage, we have the necessary ingredients to study the behavior of the linear entropy 5, 
which is rich. This is because S is actually depending on three parameters: the truncation s of the 
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WeylHeisenberg algebra, the levels of excitation \z\ of the input coherent states and the reflection 
coefficient r of the beam splitter device. We first plot the entanglement against for different finite 
dimensional Hilbert (the dimension is 2s + 1) for coherent states having an excitation level z = 5. The 
results are summarized in Figure 1: 




Figure 1: Linear entropy of coherent states as a function of the square of reflection coefficient R 

the level of excitation z = 5. 



for 



It is clear that the maximum degree of entanglement is reach when one uses a 50:50 beam splitter 
(i.e. r = l/\/2) as expected. This is independent of the dimension of the Fock space. It is remarkable 
that for r = l/\/2, the linear entropy is maximal for s = 1, which corresponds to qutrits system. 
The behavior of the entropy as function of the truncation level of the WeyHIeisenberg algebra is 
represented in Figures 2. We consider separately coherent states with weak and strong excitation level 
passing throughout a 50:50 beam splitter. 



z= 1 
z = 2 
z = 3 



Figure 2. a: Linear entropy as a function of s for weak coherent state excitation. 

In Figure 2. a, it is clearly shown that for small values of the variable 2, the linear entropy decreases 
quickly and goes to zero with increasing s. This agrees with the fact that for s large, the linear entropy 
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must vanishes as the coherent states (|2ip go to Glauber ones for ordinary harmonic oscillator, which 
does not exhibit entanglement after passing throughout a 50:50 beam splitter. In the case of stronger 
excitation, the linear entropy behaves differently as it can be seen from Figure 2.b. 



. z = 40 
. z = 30 
. z = 20 
. z= 10 



! 10 i: 14 



Figure 2.b: Linear entropy as a function of s for strong coherent state excitation. 

For instance, for z = 10, the entropy undergoes an initial increase followed by a slower decrease. It 
follows that for strong excitation levels of the coherent states, the linear entropy will reach zero for 
very large value of s. To understand the behavior of the linear entropy as function of the excitation 
level of the coherent state z, we give the Figure 3: 




Figure 3: Linear entropy as a function of coherent state excitation. 

In this Figure, we plot S as function of z for different values of s. It reflects that for weak excitation 
levels z < 5, the linear entropy increases rapidly and for higher level excitations it increases slowly. 
This explains why, for strong level excitation, the linear entropy does not go to zero for higher s 
so quickly as in the case of systems with small number of states. The separability of highly excited 
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coherent states, passing throughout a 50:50 beam sphtter, occurs for quantum systems with sufficiently 
large number of states (the dimension of the corresponding Hilbert space is very large). 

Two main features were observed through the analysis of this subsection. The first is that the 
entanglement of the output states depends heavily on \z\, more than s. For weak excitation levels z 
of the coherent states, the linear entropy decreases quickly (Figure 2. a) as s increases. This changes 
drastically for strong excitation levels and the entanglement initially increases with s (Figure 2.b). The 
second interesting feature is that for large \z\, the entanglement, after an initial increase, decreases 
slowly as we vary s. This decline is very dependent on \z\ and is lower for larger \z\. It follows 
that to see coherent states (with s large) near zero entanglement, one would need a system of higher 
excitation levels. In our numerics we considered only s such as s < 15, but we may see a significantly 
low resultant entanglement for large s but with coherent states of extremely high excitation levels. 
This reflects the resistance of the bipartite entanglement of coherent states when passed through a 
beam splitter. 

3.3 Tripartite entanglement induced by two beam splitters 

To investigate the tripartite entanglement, we now consider the action of the unitary operator 

B2 = ^2,3(^2)^1,2(^1) (34) 

on the coherent state (f2T]l . The operator B2 is obtained from ([23]) by setting k = 2. This situation 
is interesting because it allows us to study bipartite as well three particle entanglement. The input 
states are of the form jz; 0; 0) where the first mode is prepared in a coherent state \z) and the second 
and third modes are in their vacuum state. The action of B2 on the state \z, 0, 0) can be evaluated as 
follows. Indedd, using (f27|l . one verifies 



B,|n,0,0) = t ^^^^t?(.rO"-g^ (35) 
from which one obtains 

2s n n—p 

^2^,0,0) = loutput) = X^X] X] (^i'>T'^P^P)\P'P'''>^-P-p') (36) 

ri=0p=0p'=0 

where 



a(n.p,p')=MC„.''^^^mnr''^^^J^^^ (37) 

The output state is a pure state of three qudits (1, 2 and 3) system in the basis 712, n^)}- Then, 

according to [52], for the pure state |output) given by (p6|) defined on (g) (g) (dim T = 2s + 1), 
the concurrence has the form 



C(loutput)) = y [3 - {Tvpl + Trpl + Tvpl)] . (38) 

where 

Pi = Tr2,3p, P2 = Tri_3p, ps = Tri,2P (39) 
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and p = I output) (output I . This formula provides us with a measure of the tripartite entanglement of 
the coherent state (j2ip after passing through two beam splitters. We can interpret the right hand of 
the equation ()38p as follows. If we regard the pair of oscillators (23) as a single subsystem, it makes 
sense to talk about the bipartite entanglement between the first oscillator system 1 and the pair (23). 
In this respect, one can treat 1 and (23) as two systems and use the linear entropy as measure of the 
degree of entanglement between them. In the same manner, one can define the bipartite entanglement 
between the subsystems 2 and (13) and the pair 3 and (12). The linear entropies in each case are 
given by 

5i = l-Trp2, S2 = l-Txpl 53 = 1-Trpi. (40) 

It follows that the concurrence ([38]) can be expressed in terms of the sum of 5" = 81 + 82 + 8^. This 
reflects that the tripartite entanglement is related to the degree of bipartite entanglement between the 
different subcomponents of the system under consideration. 

Figures 4, 5 and 6 present the behavior of the entropies Si, 82 and 83 ([m]) in terms of the square 
of the reflection coefficients ri and r2 of the beam splitters acting in the coherent state (j2ip . We focus 
on qutrits system, i.e. s = 1 and take the level of excitation of the input coherent state as z = 5. 




Figure 4: Linear entropy 5i versus Ri = rf and R2 = r"^ for {s = 1, z = 5). 

The linear entropy 5*1 measuring the degree of bipartite entanglement between the parts 1 and (23) 
viewed as a single subsystem is plotted in Figure 4. It is important to note that 81 is maximal for 
ri = l/\/2 and r2-independent. This result can be easily shown analytically using the equations ()36p . 
([39P and (P0]l . This is not the case for the entropy 82 as it is shown by Figure 5: 
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The quantity 5*2 is associated with the linear entropy measuring the degree of bipartite entanglement 
between the subsystems 2 and (13). The entropy 5*2 is maximal for (ri = l/-v/2, "^2 = 0) or (ri = 
1, r2 = l/\/2) and vanishes for ri = or r2 = 1. Indeed, using the equation ()36p . one can simply check 
that the output state is a tensorial product of a wave function corresponding to the subsystem (13) and 
the vacuum |0) of the bosonic mode 2. Consequently, the linear entropy, associated with entanglement 
of the subsystems 2 and (13), vanishes in agreement with the results of Figure 5. Similarly, we plot in 
Figure 6 the linear entropy S3 which gives the amount of entanglement between the subsystems (12) 
and 3. 




Figure 6: Linear entropy S'3 versus Ri = rf and R2 = r"^ for {s = 1, z = 5). 

The maximal value of the linear entropy /S'3 is reached for (ri = l/-v/2, r2 = 1) or (ri = l,r2 = 1/^/2) 
and vanishes for ri = or r2 = 0. This can be easily verified using the equation (j36p . Indeed, for 
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ri = 0, the output state is separable because the action of the beam spUtters leaves invariant the state 
|2;,0,0). Similarly, for r2 = 0, the third part of the system, initially in the vacuum state, remains 
unchanged when passed throughout the beam splitters device and the output state is completely dis- 
entangled. 

At this stage, one can use the linear entropies Si, S2 and S3 to study the tripartite entanglement 
of the output state under consideration. In fact, from ()38p . it is easily seen that the square of the 
tripartite concurrence is proportional to the sum of the entropies Si, S2 and ^3. This quantity is 
plotted in Figure 7: 




Figure 7: The total Linear entropy 5 = Si + 5*2 + ^3 as functions of Ri = r\ and Ri = r| for 

{s = l,z = 5). 

We notice that for coherent states of a qutrit system (s = 1) with z = 5, the maximal tripartite 
entanglement is obtained for ri = r2 = 1 /V2. It is interesting to mention that in the particular cases 
where (ri = 0, r2 = 0), (ri = l,r2 = 0), (ri = 0, r2 = 1) and (ri = l,r2 = 1), the output states are 
given by \z, 0, 0), |0, iz, 0), \z, 0, 0) and |0, 0, —z), respectively. These states, derived from the equation 
()36p . are clearly separable. This agrees with Figure 7 from which one can see that the entropy S 
vanishes only in these particular four situations. In this subsection, we only considered qutrits. This 
study can be extended to any qudit system (s arbitrary). It should be noticed that for s large and 
high excitation levels z, no entanglement can be generated at the output. This is essentially due, as 
we discussed in the previous subsection, to the fact that in this limit all the bipartite entanglement in 
the system vanishes and hence no tripartite entanglement can be generated. Indeed, for s large the 
coherent states (pTj) reduces to Glauber ones. In this case, the action of two 50:50 beam splitters (i.e. 
in = ir2 — with a Glauber coherent states \z) incident in one port and a two particle vacuum 
|0,0) on the other ports, reduces to 

\z) |0) |0) \z/V2) \iz/2) I - z/2) (41) 
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where the notation \a) stands for the Glauber coherent states, such as 

oo 



oo „ 

n! 

n=0 



This shows that in the hmiting case s — >• oo, the resulting output state is disentangled. The results pre- 
sented in this section can be extended in different directions to investigate many features of entangled 
coherent states passing throughout a quantum network of beam splitters. 



4 Conclusion 

Using the linear entropy as measure of the bipartite entanglement, we have investigated the bipartite 
entanglement of the coherent states of a truncated Weyl-Heisenberg algebra. This algebra constitutes 
an alternative way to describe quantized electromagnetic field with finite dimensional Hilbert spaces 
instead of the usual truncated harmonic oscillator introduced by Pegg and Barnett . Subsequently, 
We have discussed the use of a quantum network of k beam splitters as experimental device to generate 
easily the coherent states associated with SU{k + 1) algebra as well as the entangled photonic states. 

We have focused on network involving one or two beam splitters. In the first case, this device 
is used to investigate the degree of entanglement of coherent state associated with the truncated 
harmonic oscillator. We have shown that, for weak level of excitation, the linear entropy goes faster 
to zero for s increasing. For stronger excitation levels, the entropy decreases slowly as the Fock space 
dimension increases. The second case involves two beams splitters. In this case taking advantage 
from the relation of the tripartite concurrence and the partial linear entropies (see equation (j38p ). we 
have discussed the bipartite as well tripartite entanglement of a coherent states passing throughout a 
network of two beam splitters. Moreover, we have shown that the maximal tripartite entanglement is 
reached when the beam splitters are 50:50. 
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